Finite-size scaling as a tool in the search for the QCD critical point in heavy ion data 
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Given the short hfetime and the reduced volume of the quark-gluon plasma (QGP) formed in 
high-energy heavy ion collisions, a possible critical endpoint (CEP) will be blurred in a region 
and the effects from criticality severely smoothened. Nevertheless, the non-monotonic behavior of 
correlation functions near criticality for systems of different sizes, given by different centralities in 
heavy ion collisions, must obey finite-size scaling. We apply the predicting power of scaling plots to 
the search for the CEP of strong interactions in heavy ion collisions using data from RHIC and SPS. 
The results of our data analysis exclude a critical point below chemical potentials /i ~ 450 MeV. 
Extrapolating the analysis, we speculate that criticality could appear slightly above ~ 500 MeV. 
Using available data we extrapolate our scaling curves to predict the behavior of new data at lower 
center-of-mass energy, currently being investigated in the Beam Energy Scan program at RHIC. If 
it turns out that the QGP phase is no longer achievable in heavy ion experiments before the CEP is 
reached, FSS might be the only way to experimentally estimate its position in the phase diagram. 
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High-energy heavy ion collisions can produce energy 
densities that are high enough to generate a new state of 
matter, the quark-gluon plasma (QGP), probing a new 
region in the phase diagram of strong interactions and 
testing quantum chromodynamics (QCD) under very ex- 
treme conditions. These experimental events are often 
referred to as little bangs since they reproduce some of 
the conditions found in the early universe at the time 
of the primordial quark-hadron transition. Nevertheless, 
one should bear in mind that the space-time scales that 
set the dynamics of the quark-gluon plasma formed in 
heavy ion collisions differ from the cosmological ones by 
almost twenty orders of magnitude. This means that the 
early universe at the time of the QCD phase transition(s) 
was already very large, in the sense that a description as- 
suming the thermodynamic limit is well justified, whereas 
in the case of heavy ion collisions one should be more 
cautious. 

In fact, even if this issue is usually overlooked, it has 
been shown using lattice simulations [l|, 0] and different 
effective model approaches that finite-size effects 

can dramatically modify the phase structure of strong in- 
teraction, dislocating critical lines and critical points, and 
also affect significantly the dynamics of phase conversion 
0, Q. Moreover, as discussed in Ref. p], the modifica- 
tions suffered by the phase diagram features, such as the 
critical line and isentropic trajectories, are sensitive to 
the boundary conditions, as expected for small systems. 
Therefore, the data corresponding to small QGP systems, 
such as the ones formed in non-central heavy ion colli- 
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sions, should not be naively compared to what one would 
expect in the thermodynamic limit. A direct consequence 
of this fact is that all signatures of the second-order criti- 
cal endpoint based on the non-monotonic behavior 
or sign modifications ll| of particle correlation fluctua- 
tions will probe a pseudocritical endpoint that can be 
significantly shifted from the genuine (unique) critical 
endpoint by finite-size corrections and will be sensitive 
to boundary effects. This feature, together with the even 
more crucial limitation on the growth of the correlation 
length due to the finite (short) lifetime of the plasma 
state and critical slowing down [l^, makes the ex- 
perimental searches of signatures of the presence of a 
critical point at lower energies very challenging, even if 
one considers higher moments which grow faster with the 
correlation length as suggested in Ref. The round- 

ing and smoothening of fluctuation peaks tend to hide 
them behind the background. Besides, measuring higher 
moments of fluctuations of experimental observables in 
heavy-ion collisions is also a challenging task. 

On the other hand, given the fact that the thermal en- 
vironment corresponding to the region of QGP formed 
in heavy ion collisions can be classified according to the 
centrality of the collision, so that events can be separated 
according to the size of the plasma that is created, the 
finiteness of the system also brings a bright side: the pos- 
sibility of finite-size scahng (FSS) analysis (For 
reviews, see Refs. [iiSl). FSS is a powerful statistical 
mechanics technique that prescinds from the knowledge 
of the details of a given system; instead, it provides infor- 
mation about its criticality based solely on very general 
characteristics. Although it is clearly not simple to de- 
fine an appropriate scaling variable in the case of heavy 
ion collisions, the flexibility of the FSS method allows 
for a pragmatic approach for the use of scaling plots in 
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the search for the critical endpoint as was dehneated in 
Ref. The essential point is that although the re- 

duced volume of the QGP formed in high-energy heavy 
ion collisions will dissolve a possible critical point into 
a region and make the effects from criticality severely 
smoothened, the non-monotonic behavior of correlation 
functions for systems of different sizes, given by different 
centralities, must obey FSS near criticality 19l421|. 

In this paper we apply the predicting power of scal- 
ing plots to the search for the CEP of strong inter- 
actions in heavy ion collisions using data from RHIC 
and SPS. Defining appropriate scaling variables, we 
generate scaling plots for data sets with ^/snn = 
17.3,19.6,62.4,130,200 GeV. From the results of our 
data analysis we can exclude a critical point below chem- 
ical potentials fi ~ 450 MeV. This seems to be consistent 
with lattice results that exclude a chiral critical point in 
the region of fi smaller or comparable to the temperature 

Extrapolating the analysis, we speculate that critical- 
ity could appear slightly above n ^ 500 MeV. Using avail- 
able data we extrapolate our scaling curves to predict 
the behavior of new data at lower center-of-mass energy, 
currently being investigated in the Beam Energy Scan 
program at RHIC. If it turns out that the QGP phase is 
no longer achievable in heavy ion experiments before the 
CEP is reached, FSS might be the only way to experi- 
mentally estimate its position in the phase diagram. 

Finding the appropriate scaling variable(s) is almost 
never an easy exercise for real physical systems. In this 
analysis we pragmatically choose to work with a one- 
dimensional scaling function motivated by the simplifi- 
cation of thermal model descriptions of the freeze-out re- 
gion [l8| . which connects temperature and chemical po- 
tential, so that our scaling variable can be built from 

It is well known that as one goes near criticality in 
a second-order phase transition the correlation functions 
of the system scale with some power of the correlation 
length. This feature is independent of the details of the 
system under consideration, and the power is dictated by 
which correlation function one considers and the univer- 
sality class to which the system belongs 1^, 23|. There- 
fore, it is clear that observables that can be written in 
terms of the correlation functions will exhibit a non- 
monotonic behavior near a critical point, the effect being 
more dramatic as one goes up in moments of the fluc- 
tuations. This behavior provides a clean evidence of the 
presence of a critical point Thus, it could be used 

as a precise signature in the analysis of data from heavy 
ion collisions using event- by- event analysis, as proposed 
in Refs. were it not the case that these systems are 

usually very small and come in different sizes. Therefore, 
the peaks that are typical from this divergent behavior 
in the thermodynamic limit will be smoothened by the 
effects of finite size and short lifetime of the systems un- 



der consideration in heavy ion collisions. Even higher 
moments which could provide stronger divergences 
(that grow with higher powers of the correlation length) 
are more difficult to analyze under these circumstances, 
from the experimental point of view. So, as an alternative 
in data analysis, besides looking for non-monotonicity in 
particle multiplicities, one should also test for the scaling 
that differently-sized systems must obey near criticality, 
i.e. finite-size scaling. This also corresponds to the cor- 
rect way to extrapolate the results of measurements of 
the finite system to the thermodynamic limit. 

The FSS hypothesis [l^ was conjectured before the de- 
velopment of the renormalization group (RG). However, 
it can be derived quite naturally by applying RG tech- 
niques to critical phenomena [16j, |21j . For a system with 
typical linear dimension L, the singular part of the free 
energy density scales as 



M{g},L'^)^£-'^fs{{9'},iL-^) 



(1) 



when the lengths are reduced by a factor £ in a given 
RG transformation, {g} are couplings and d is the space 
dimension. The fact that the system is finite is irrele- 
vant for this implementation, since the RG transforma- 
tions are local. Near a RG fixed point, one can write ([T]) 
in terms of the right eigenvectors of the linearized RG 
transformation: 



fs{t,r,...,L-^)^r'ifs{t£y\rey^...,£L-^) 



(2) 



so that one notices that L^^ behaves like a relevant eigen- 
vector with eigenvalue A^, = £, and — 1. Here, 
t = (T — Tc)/Tc is the reduced temperature, T^. being 
the temperature associated with the critical point. The 
reduced temperature f is a dimensionless measure of the 
distance to the critical point when no other external pa- 
rameter is considered. In the presence of other external 
parameters, such as the baryonic chemical potential /z, 
one should redefine this distance accordingly to include 
the dependence on r = (/i — /ic)//ic7 which plays a role 
analogous to the magnetic field variable in a Ising system. 
Tc and /ic are defined in the thermodynamic limit, and 
the true second-order phase transition occurs for vanish- 
ing t, r and (and other couplings reaching their value 
at the fixed point, {g = <?*}). 

For finite L, crossover effects become important. If 
the correlation length diverges as ^oo ~ t~'^ at criticality, 
where v is the corresponding critical exponent, in the 
case of L^^t^'^ ^ 1 the system is no longer governed by 
the critical fixed point and L limits the growth of the 
correlation length, rounding all singularities [26|. If L 
is finite, ^ is analytic in the limit t — >■ 0, etc and one 
can draw scaling plots of L/^ vs. a given coupling g for 
different values of L to find that all curves cross at g = g* 
in this limit, which is a way to determine g* . The critical 
temperature can also be determined in this fashion, since 
the curves will also cross at t = tn- 
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This scahng plot technique can be extended, taken to 
its full power for other quantities, such as correlation 
functions. An observable X in a finite thermal system 
can be written, in the neighborhood of criticality, in the 
following form [l6]: 



X{t,{g}-l-L) = L'<^'^f(tL^/'^) , 



(3) 



where is the bulk (dimension) exponent oi X and {g} 
dimensionless coupling constants. The function /(?/) is 
universal up to scale fixing, and the critical exponents are 
sensitive essentially to dimensionality and internal sym- 
metry, which will give rise to the different universality 
classes 3, 21 1 . To simplify the discussion, we ignore for 
the moment the chemical potential variable r. In prin- 
ciple, / is a function of two scaling variables, though. 
Again, if one plots the scaled observable versus i, the 
curves should cross aX t = tc- Moreover, using the ap- 
propriate scaling variable, instead of all curves should 
collapse into one single curve if one is not far from the 
critical point. So, this technique can be applied to the 
analysis of observables that are directly related to the 
correlation function of the order parameter of the transi- 
tion, such as fluctuations of the multiplicity of soft pions 
i- 

This is the idea of using full scaling plots, a method 
that has been proven to work even for tiny systems, to 
search for the critical endpoint of QCD in heavy ion 
data. Here, as in statistical mechanics searches for criti- 
cal points, one can treat Tc, /ic and the critical exponents 
as parameters in a scaling plot fit, the system size being 
provided by different centrality bins or by the number 
of participants in the collision. In fact, the scaling vari- 
able is defined up to L- and ^-independent multiplicative 
factors, so that the knowledge of the actual size of the 
system is not needed. 

As mentioned previously, the correct scaling variable 
should measure the distance from the critical point, 
thereby involving both temperature and chemical poten- 
tial, or their reduced versions t and r. This produces a 
two-dimensional scaling function and makes the analy- 
sis of heavy ion data highly nontrivial. Phenomenologi- 
cally, we adopt a simplification motivated by results from 
thermal models for the freeze-out region, connecting tem- 
perature and chemical potential. We can parametrize the 
freeze-out curve by ^/sWn, and build our one-dimensional 
scaling variable from this quantity and the size of the sys- 
tem. 

The range of sizes that can be accessed in heavy ion 
collisions over which FSS can be tested is limited. It 
is reasonable to assume that for a locally thermalized 
quark gluon plasma to form, the system needs to be sev- 
eral times the hadronic size (several fm). The largest 
possible system will have a diameter of approximately 
15 fm. Estimates from an analysis of HBT data from 
STAR indicate that from peripheral to central collisions, 
the system size changes by a factor of 3 — 4 22| . We plot 



the scaled observable vs [{fi — fic)/ Hc] L^^'^ ■ The value 
of L for the various centralities can be estimated from 
a Glauber Monte-Carlo model. For our study, the over- 
lap area is calculated from S± = ~ sinO) where 
Q ~ 2cos''^{b/2RAu), where Rau is the nuclear radius 
for Au and b is the collision impact parameter. Then the 
length is taken to be L = 2yJS_\_/'K. 

A system of size Li can be compared to a system of 
size L2 only when 



^1 - = (a*2 ^ A*c) ( ^ 



(4) 



This constrains which /i values can be compared when 
testing for scaling by directly comparing data rather than 
by extrapolations. Since a central Au + Au collision is 
about four times larger than a peripheral one {L2 — 4Li) 
and taking v = 2/3, we find that a measurement in pe- 
ripheral Au-l-Au collisions at /ii can be compared to a 
central collision at ^2 when /ii — /i^ ~ 8(/i2 — /ic)- 

Data on various fluctuations have been measured at 
RHIC and the SPS over an energy range from ^^snn = 5 
GeV to 200 GeV and the centrality dependence of pr fluc- 
tuations have been measured at RHIC for energies from 
19.6 GeV to 200 GeV. This data provides an opportunity 
to look for evidence of finite-size scaling. 

To search for scaling, we consider the correlation mea- 
sure (Tpy /(pt) [23 scaled by L-'i-/", according to Eq. H 
We consider the pt fluctuations CTp^ scaled by {pt) to 
obtain a dimensionless variable ^ . We use the correlation 
data measured in bins corresponding to the 0-5%, 5-10%, 
10-20%, 20-30%, 30-40%, 40-50%, 50-60%, and 60-70% 
most central collisions. We estimate the corresponding 
lengths L to be 12.4, 11.1, 9.6, 8.0, 6.8, 5.6, 4.5, and 3.4 
fm. The exponent v — 2/3 is determined by the Ising 
universality class of QCD and we consider values of 
around 1 (ignoring small anomalous dimensions correc- 
tions). We also varied the value of from 0.5 to 2.0 
and found that changing ^x within this range does not 
improve the scaling behavior. Using a smaller value of ^x 
simply causes Gp^/ipT) scaled by L~^^/^ to drop more 
slowly with increasing centrality and increasing ^x causes 
this drop to happen more quickly. 

In Figs.[T]and[2]we plot ap^/{pT) scaled by L^f^/'^ vs fi 
for different system sizes, using data extracted from col- 
lisions at ^/sWn = 19.6, 62.4, 130, 200 GeV. We use the 
parameterization from Ref. |25j to convert from ^/snn to 
fi. If there is a critical point at /x = fj-cnt, the curves for 
different sizes of the system should cross at this value of ^. 
However, since the currently available data is restricted 



^ The STAR collaboration has also published pT fluctuations us- 
ing an alternative variable Acp^, |24( . The results found us- 
ing this variable arc qualitatively equivalent to the ones for 
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FIG. 1: Scaled o-pt/{pt) vs jj. for diflerent system sizes, and 
with 1/ = 2/3 and 72, = 1. Data extracted from RHIC colli- 
sions at energies ^/snn = 19.6, 62.4, 130, and 200 GeV (linear 
fit, see text). 

to not so large values of the chemical potential, one has to 
perform extrapolations using fits. The scaling function / 
in Eq. [3] is expected to be smoothly varying around the 
critical point, so we fit the data corresponding to a given 
linear size L to a polynomial, but constraining the poly- 
nomials to enforce the condition that all the curves cross 
at some /i = /icrit, where /iciit is an adjustable parameter 
in the fit. This clearly assumes the existence of a critical 
point. In Fig. [1] we use a linear fit. The approximate en- 
ergy independence of (jp^/{pT) along with the linear fit, 
leads to a very large /i value where the curves can cross 
(/Lt ~ 3 GeV). There's no reason however to assume a 
linear fit function, so in Fig. [5] we also try a second order 
polynomial. Using the a second order polynomial func- 
tion for / allows the curves from different system sizes 
to cross at a much smaller value of /i. Based on this fit, 
we find that the data is consistent with a critical point 
at /X ~ 510 MeV corresponding to a ^smn of 5.75 GeV. 

The value estimated for /^crit based on finite-size scal- 
ing of current data is highly dependent on the assumed 
functional form of /. The approximate energy indepen- 
dence of (Tpy/(pT) for a given L, however, already indi- 
cates within the finite-size scaling assumption that the 
critical point should be at values well above those cur- 
rently available. Based on finite-size scaling of apj./{pT), 
one would not expect a critical point at /z < 400 MeV. 
Having data for lower values of i/s, i.e. higher values of 
chemical potential, as expected from the analysis of the 
Beam Energy Scan program at RHIC [23|, one should 
be able to study full scaling plots of ap^/{pT) scaled by 
l^-i^/v p--Pcrit j^i/i/ .^vithout the need of long extrap- 
olations. For the current set of data, these full scaling 
plots are still not very enlightening. 

RHIC has also run at lower energies in order to search 
for a critical point in the Beam Energy Scan program. 
That data is currently being analyzed. Here we use the 
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FIG. 2: Scaled ap^/{jpT) vs /i for different system sizes, 
assuming /icrit = 509 MeV, which corresponds to a crit- 
ical point at y'sjvjv ~ 5.75 GeV. Again, v — 2/3 and 
73; = 1. Data extracted from RHIG collisions at energies 
y/SNN ~ 19.6, 62.4, 130, and 200 GeV (second order polyno- 
mial fit, see text). 

quadratic polynomial fit of STAR data (Fig. [2]) and as- 
sume the critical point is at 509 MeV to make predictions 
for apj./{pT) at lower energies. The finite-size scaling 
scenario along with currently available data, allows us to 
predict the energy and system size dependence of fluctu- 
ations for any given values of /^crit: Ix and v. We show 
this expectation as a function of the number of partici- 
pants, A'^pait, for three proposed beam energies: 11.5,7.7 
and 5 GeV in Fig. |31 Notice that the centrality depen- 
dence changes once one moves to the other side of the 
critical point. This is a condition enforced by finite-size 
scaling which provides a generic signal for having reached 
the first-order phase transition side of the critical point. 
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FIG. 3: The expected measurement of (Tp^/ {pt) as a function 
of the number of participants at lower energies assuming the 
critical point is at 509 MeV as extracted from the quadratic 
polynomial fit of STAR data. 

Searching for the the critical endpoint in the phase di- 



5 



agram for strong interactions is a remarkably challenging 
task. On one hand, lattice simulations can not provide 
the sort of guidance that is possible at zero density, where 
there is no sign problem, and effective models point to its 
existence, qualitatively, but yield very different quantita- 
tive predictions. On the other hand, heavy ion collision 
experiments probe limited regions of the phase diagram. 
Besides, they come with non-trivial background contri- 
butions that tend to blur the signatures provided by the 
non-monotonic behavior of observables built from corre- 
lation functions of the order parameter. 

The fact that finite-size scaling prescinds from the 
knowledge of the details of the system under consider- 
ation, providing information about its criticality based 
solely on its most general features, makes it a very pow- 
erful and pragmatic tool for data analysis in the heavy ion 
collision search for the critical point. From a very limited 
data set in energy spam, we have used FSS to exclude the 
presence of a critical point at small values of the baryonic 
chemical potential, below 450 MeV. We have also used 
the scaling function to predict the behavior of data with 
system size at lower energies. We are looking forward to 
compare our predictions to the outcome of data analysis 
from the Beam Energy Scan program at RHIC. 
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